simulation results/work_for_gamma_results.R

In BriceonWiley/IntegratedLikelihood.R:

pics_coverage <-
  list(
    results_050, results_100,
    results_150, results_200,
    results_250
  ) %>%
  map(~{
    map_dfr(.x, ~ coverage_gam(.x, TRUE)) %>%
      select(-N, -NaNs) %>%
      gather("Interval", "Coverage", -gamma, -p1, -p2) %>%
      ggplot(aes(gamma, Coverage)) +
      geom_path(aes(linetype = Interval), size = 1) +
      geom_hline(yintercept = 0.95, linetype = "dashed") +
      theme_minimal(25) +
      ylim(c(0.6, 1)) +
      scale_x_continuous(name = expression(gamma))
  })

pics_width <-
  list(
    results_050, results_100,
    results_150, results_200,
    results_250
  ) %>%
  map(~{
    map_dfr(.x, ~ width_gam(.x, TRUE)) %>%
      select(-N, -NaNs) %>%
      gather("Interval", "Width", -gamma, -p1, -p2) %>%
      ggplot(aes(gamma, Width)) +
      geom_path(aes(linetype = Interval), size = 1) +
      theme_minimal(25) +
      ylim(c(0, 6)) +
      scale_x_continuous(name = expression(gamma))
  })

map2(
  pics_coverage, str_pad(seq(50, 250, 50), 3, "left", 0),
  ~ ggsave(glue::glue("c03_coverage_{.y}.pdf"), .x)
  )
map2(
  pics_width, str_pad(seq(50, 250, 50), 3, "left", 0),
  ~ ggsave(glue::glue("c03_width_{.y}.pdf"), .x)
  )